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We explore how the local environment is related to the redshift, type, and 
luminosity of active galactic nuclei (AGN). Recent simulations and observations 
are converging on the view that the extreme luminosity of quasars is fueled in 
major mergers of gas-rich galaxies. In such a picture, quasars are expected to be 
located in regions with a higher density of galaxies on small scales where mergers 
are more likely to take place. However, in this picture, the activity observed in 
low-luminosity AGN is due to secular processes that are less dependent on the 
local galaxy density. To test this hypothesis, we compare the local photomet- 
ric galaxy density on kiloparsec scales around spectroscopic Type I and Type II 
quasars to the local density around lower luminosity spectroscopic Type I and 
Type II AGN. To minimize projection effects and evolution in the photometric 
galaxy sample we use to characterize AGN environments, we place our random 
control sample at the same redshift as our AGN and impose a narrow redshift 
window around both the AGN and control targets. We find that higher luminos- 
ity AGN have more overdense environments compared to lower luminosity AGN 
on all scales out to our 2 hy Q Mpc limit. Additionally, in the range 0.3 ^ z ^ 0.6, 
Type II quasars have similarly overdense environments to those of bright Type I 
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quasars on all scales out to our 2 hj Mpc limit, while the environment of dim- 
mer Type I quasars appears to be less overdense than the environment of Type 
II quasars. We see increased overdensity for Type II AGN compared to Type 

1 AGN on scales out to our limit of 2 Mpc in overlapping redshift ranges. 
We also detect marginal evidence for evolution in the number of galaxies within 

2 /ifo 1 Mpc of a quasar with redshift. 

Subject headings: cosmology: observations — large-scale structure of universe - 
quasars: general — surveys 



INTRODUCTION 



The local environments of quasars provide y aluable insight into the formation history 



an d evolution of AGN (e.g., 



E 



lin ggon et al.lll99ll ). Quasar environments were first studied 



by Bahcall. Schmidt, fe GunnT ( 1969), who used a sampl e of five quasars to show that quasars 
are associated with galaxy clusters. lYee fc Greenl (11984 ) found that quasars reside in regions 



with higher galaxy density, and more recent work has confirmed that quasars are found in 



regions of galaxy groups or clusters of poor to moderate richness (jBahcall fc Chokshi 



Fisher et al.lll996l : McLure fc Dunlopll200ll : IWold et al.ll200ll : IColdwell et al 



2002 



1991 



Barr et al. 



20031 ) . Although studies of several X-ray- and radio-selected s amples have found e vidence for 



a relationship betwe en environment and AGN activity (e.g., iWurtz et al.l 119971 ; iBestl 12004 ; 
Sochting et al.ll2004l ). the Sloan Digital Sky Survey (SDSS) is the first survey to allow mean- 
ingful studies of quasar environments, be cause it samp l es lar ge numbers of both quasars and 
galaxies at z < 0.4. Using SDSS data, ISerber et al.l (120061 ) concluded that the density of 
photometric galaxies around quasars increases with decreasing angular scale, but is indepen- 
dent of redshift for z ^ 0.4. They also provided evidence for a higher density of galaxies 
around more luminous quasars at scales less than 100 Yr x kpc, wh ile at larger angular scales 



the density appears to be largely indepe ndent of luminosity (c.f. iPorciani fc Norbergi 12006 
da Angela etaDhoQsl : levers et aliboOTT ). 



The results of ISerber et alj (2006) agree with other studies showing enhanced clustering 



of quasars on small scales. iDjorgovskil (Il99ll ) first linked the excess of quasar clustering on 
small scales to galaxy interactions. Studies of small-scale clustering of quasars (e.g., binary 
and triplet quasars) also support the hypothes i s that there is excess quasar clustering on 



scales of < 100 h 1 kpc ( 


Kochane 


<l et al. 


1999 


Diorgovski et al. 


2007; 


Mvers et al. 


2007|). 



1999; Mortlock et al. 1999; Hennawi et al. 2006 



An excess of quasar pairs on small scales naturally follows from a merger origin for quasar 
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activity, whether the se pairs simply trace biased groups whe re mergers are likely to occur 



(Ho pkins et al.ll2008f) or are being excited in merging galaxies (lDjorgovskilll99ll ; iMyers et al. 



20081 ). iHopkins et al.l (120061 ) have developed a unified, merger- driven framework th at natu- 
rally predicts that quasar environments should be highly biased (IHopkins et al.ll2008l ). These 
simulations show that major mergers between gas-rich galaxies are the likely mechanisms to 
trigger bright quas ar activity, and that this a ctivity is a phase in the evolution of massive 
spheroidal galaxies (IHopkins et al.ll2005l . 120081 ). In contrast, secular mechanisms may fuel the 
activity in most low-luminosity A GN, implying that the small-scale environments of these 
objects should have a smaller bias (IHopkins fc Hernquistll2006l ; IHopkins et al.ll2008l ). There- 
fore, we would expect that objects driven by major mergers will have biased environments 
on small scales, whereas objects fueled by secular means will reside in less rich environ- 
ments. Such a simplification hides many subtleties, however, as secular mechanisms such as 
harassment can probably only occur in slightly overdense environments. Further, for objects 
whose observed c haracterist ics differ purely because of viewing angle or internal structure 
( )Antonuccilll993l ; lElvid |2000| ) . there should be no particular difference in local environment. 
This, of course, would only be the case if that structure is not correlated with fueling, as 
could occur, for instance, if more luminous quasars had strong winds. It is therefore impor- 
tant to understand the relationship between the physical properties of AGN and their local 
environment, which will in turn provide insight into what aspects of AGN properties are 
explained by formation history, fueling, or simply by structure and orientation. 

In this paper, we address this challenge b y studying th e natu re of AGN environments. 
We improve upon the most recent SDSS study, ISerber et al.l (120061 ). in several ways. We use 
larger samples of background photometric galaxies, as well as larger, more recent samples of 
spectroscopic AGN. Our spectroscopic data is divided into four target samples: Type I and 
Type II quasars (e.g., AGN with the highest intrinsic luminosity) and lower- luminosity Type 
I and Type II AGN. The spectra of Type I AGN and quasars are characterized by broad 



emission lines (FWHM > 1000 km s -1 ; e.g., lHao et al.ll2005al : ISchneider et all 120071) while 



the spectra of Type II AGN and quasars exhibit narrow emission lines (e.g.. lHao et al.ll2005al ; 
Zakamska et al.ll2003l ). We compare the overdensity of Type I quasars to the overdensity 
of Type II quasars as well as to lower luminosity Type I and Type II AGN. Additionally, 
we include cuts in photometric redshift space around spectroscopic targets and the random 
positions to which they are compared to minimize interloping foreground or background 
objects, as well as marginalize over any redshift evolution of the photometric galaxy sample. 
By using photometric redshift cuts, we obtain more realistic overdensity estimates and errors, 
and we are able to extend the study of AGN environments in the SDSS to z ^ 0.6. 

The data samples used in this paper are discussed in Section [21 and the technique 
for calculating overdensities is described in Section [31 We give a detailed discussion of the 
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relationship of type, redshift and luminosity with the local environment of AGN in Section [U 
Throughout the paper, we assume a concordance cosmology Om = 0.3, Qa = 0.7, H = 70 
km s _1 Mpc -1 (with h = 0.7) in order to compare to results from previous studies. 



2. DATA 

We study the environments of spectroscopic AGN targets by counting photometric 
galaxies within a 2.0 hj Mpc projected comoving distance of the target center (i.e., we 
consider a conical slice around the target rather than a spherical volume). The samples 
of Type I quasars, Type II quasars, and lower luminosity AGN, and photometric galaxies 
selected from the SDSS are described below. We apply masks to the spectroscopic samples 
to eliminate objects within 2.0 Mpc of the survey edge or an area masked out by SDSS 
or with seeing > 1.5" or r-band reddening A r ^ 0.2. Figure [1] shows the redshift distribution 
of the spectroscopic and photometric samples used in our analysis. 



2.1. Spectroscopic Targets 



Our sample of spectroscopic Type I quasar tar gets is drawn fro m the the SPSS Fift h 
Data Release (DR5; lAdelman-McCarthy et al. I l2007h Quasar Catalog (ISchneider et al.ll2007l ) , 
which includes i^-corrected absolute i-band magnitudes for each object. Our quasar sample 
has —27.5 ^ Mj ^ —22.0, and our masking reduces the sample to 6, 034 quasar targets with 
0.11 ^ z ^ 0.6. In addition to absolute magnitude information, the catalog provides a flag to 
distinguish between resolved and point source objects. The measured luminosity of resolved 
sources will likely be contaminated by starlight from the host galaxy, so we have excluded 
the 1, 800 quasars with extended morphology. The final Type I quasar sample contains 4, 234 
objects. 



We draw Type II quasar targets from the sample presented by IZakamska et al.l (120031 ). 
After cutting the sample to match the high redshift limit of the main quasar sample [z ^ 0.6) 
and masking this sample as described above, we have 160 targets with 0.3 ^ z ^ 0.6. 

We compare the quasar targets to lower-luminosity Type I and Type II AGN from Hao 
(private communic ation) sele c ted fro m SDSS DR5 spectroscopic galaxies according to the 
criteria laid out in lHao et al.l (j2005al ). The class ification of these galaxies as AGN depend s 
on the strengths of the [O III] and B./3 lines ( lHao et al.l l2005al ; IKauffmann et al.l 120031 ) ; 
therefore our low-redshift limit is set to z = 0.1 1, as this is where the [O III] (AA4959, 5007) 



lines enter the r-band (IKauffmann et al.l 120031 ) . resulting in a more uniform classification. 
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After mas king as above, there are 1, 744 Type I AGN and 3, 407 Type II AGN following the 
criteria of iKewley et al.l (120011 ) . While we could have adopted the less stringent criteria of 



Kauffmann et al.l (120031 ). we wish to be conservative in our s ample selection in this analysis 



and minimize the contribution of non-accretion luminosity (lHao et al.ll2005al ). The lower- 
luminosity AGN samples have a redshift range of 0.11 ^ z ^ 0.33, with the majority of 
sources at z < 0.15. 



2.2. Photometric Galaxies 



The photometric galaxies are drawn from the SDSS DR5 database by selecting all pri- 
mary objects photometrically classified as galaxies with r-band extinction corrected magni- 
tude in the range 14.0 < r < 21.(0 All of these obje cts have been assign ed photometric 
redshifts via a template- fitting technique as described in lCsabai et al.l (120031 ) . Our final pho- 
tometric galaxy sample of over 28 million (28, 856, 324) objects consists of only those objects 
that pass the flag requirements for a clean galaxy sample! While we do not make explicit 
redshift cuts on the photometric galaxy sample, our technique (as described in Section 13.11) 
effectively limits the sample to 0.06 ^ z ^ 0.65. 



3. TECHNIQUE 

We count the number of photometric galaxies within a comoving radius of 2.0 Mpc 
of each spectroscopic target (e.g., spectroscopic quasar, AGN, or spectroscopic galaxy), ex- 
cluding any galaxies that are within 25 h^Q kpc of the target. At z < 0.4, 25 h^Q kpc 
correspo nds to an angular size of > 3.3' ', which is approximately twice the average seeing 
in DR5 (lAdelman-McCarthy et al.ll2007l ). At angular scales smaller than this, deblending 
complicates the reliable detection and measurement of faint galaxies. 

We generate a large number of random positions in the DR5 footprint area for each 
redshift increment of 0.001 in our redshift range of 0.11 ^ z $C 0.6. We mask the positions 
in the same manner as we masked the spectroscopic targets, leaving at least 1,000 random 
positions that are more than 2.0 Mpc away from the survey edge or a masked area for 
each redshift value. We note that this approach helps counteract any bias due to evolution in 



1 Note that Serber et al. ( 2006h used an r-band limit on their photometric galaxies rather than i-band as 
stated in their paper (W. Serber and R. Scranton, private communication). 



as defined by http://cas.sdss.Org/astro/en/help/docs/realquery.asp#flags 
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the photometric galaxy sample (see also Section I37L]) . We count the number of photometric 
galaxies within a designated comoving distance around random positions and calculate the 
mean cumulative number of counts for that redshift increment as 



Rj 



N z 



where N z is the number of random postions R z at a given redshift increment z. We calculate 
the error corresponding to the mean random counts as 



rFfi tt 



a 



Ri 



(2) 



which is the variance on the mean random counts at a given redshift increment z. 



The cumulative bincounts Cj around spectroscopic targets are matched with the mean 
cumulative random bincounts Ri (and error e^J at the redshift increment closest to the 
target's redshift. We calculate a mean overdensity 5bm in a particular scale, redshift, or 
absolute magnitude bin as 



'bin 



N 



- 1 



a 



bin 



N 



Rhi 



- 1 



Chi 



bin 



bin 



R 



- 1 



(3) 



bin 



where Cj is the counts around each target in the bin, Ri is the mean counts around random 
positions at the corresponding redshift, and there are TV total targets in the bin. The error 
on the overdensity is determined via error propagation: 



p 2 r" 2 

D2 D4 C Rbin 

^bin n ^ 



(4) 



bin 



where e Cbm = ^fQ~ and el = ^ e 



N „2 
Ri 



We will refer to the quantity of as t he mean density; this quantity is used to compare 



our results to those of Serber et al. 



(120061 ). In order to compare the environments of our 



various spectroscopic populations, we use the mean overdensity, — 1, which is associated 
with the underlying dark m atter distribution and can be more directly related to correlation 
analyses of clustering (e.g.. IPadmanabhan et al.ll2008l ). 



3.1. Sz Cut 



One of the difficulties in using photometric galaxy samples for overdensity measurements 
is the issue of projection effects, where foreground or background objects contaminate a 
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measurement. We use photometric redshifts assigned to the photometric galaxies to minimize 
this complication. We apply a photometric redshift cut on the galaxies so that only those 
galaxies which satisfy \z tar get — z p hotogai\ ^ 8z are counted in each bin. Crucially, the same 
5z cut is applied to both the spectroscopic targets and the random positions to which they are 
compared, which, as noted above, are also placed at the same redshift as the spectroscopic 
targets. We therefore minimize contamination by most galaxies outside of the Sz interval. 
Additionally, by calculating our spectroscopic-photometric and random-photometric counts 
in the same z ± 5z bin, we marginalize redshift evolution in the photometric galaxy sample 
outside of that z ± 5z bin. We assume for the purposes of this work that there is no 
redshift evolution in the photometric galaxy sample over this small 5z interval. We have not 
accounted for the changes in photometric redshift accuracy as a function of magnitude and 
redshift, which we assume are negligible for this analysis. 

We verify that the projection effect issue is mitigated by the 5z cut without introducing 
systematics by calculating overdensities for random positions with the same redshift distri- 
bution as the Type I quasar sample. We find that the overdensities of photometric galaxies 
around random positions is consistent with zero on all scales with and without the photo- 
metric redshift cut. We use the value 8z = 0.05, which is large enough to encompass the 



effect ive rms error of the photometric redshifts (Az rms = 0.04 for r < 18; iBudavari et al. 



20031 ). for all further analysis. We have tested other values of 5z and find that they give 
consistent results, albeit with larger uncertainties for narrower cuts, which is consistent with 
the expectations of Poissonian sampling. 



To compare directly to the results of ISerber et al.l (120061 ). who did not apply any such 



redshift cut, we calculate the mean density of photometric galaxies around Type I quasars 
with -24.2 ^ M t ^ -22.0 and 0.08 < z < 0.4. At a scale of 250 kpc, the density of 
photometric galaxies around quasars is 1.41 ± 0.033 and around L* galaxies is 1.15 ± 0.005 
without the 5z cut. However, applying the Sz decreases the random background noise, so 
we measure an environment density of 2.11 ± 0.096 around quasars and 1.74 ± 0.020 around 
L* galaxies at the same scale. In order to confirm that we have not added systematics by 
using the 5z cut, we compare the relative densities of quasar environments to L* galaxy 
environments. The relative density of photometric galaxies around quasars to that around 
L* galaxies is 1.22 ± 0.029 without the Sz cut. The relative density does not appreciably 
change when the 5z cut is used: we find the relative density to be 1.22 ± 0.057. The true 
physical effect of the 5z cut is shown in the comparison of mean overdensities. At the 
same scale of 250 /ifg 1 kpc, the relative overdensity around quasars compared to around L* 
galaxies is 2.67 ± 0.236 without the Sz cut, but is 1.51 ± 0.137 with the Sz cut. Because we 
have removed projection effects, the relative overdensities are lower when the Sz cut is used; 
however, the errors on the mean densities with the 5z cut have increased. We believe these 
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larger errors are more physically relevant: with no 5z cut, objects not actually correlated with 
the target will reduce Poissonian error estimates. We use the Sz cut to extend our redshift 
range to include spectroscopic targets with 0.11 ^ z ^ 0.6 without concern that foreground 
objects will contaminate the overdensity measurements. Therefore, all subsequent analysis 
and figures include the Sz = 0.05 cut. 



3.2. Measurement 

In Figure EJ we present the mean cumulative overdensity of photometric galaxies as a 
function of scale for each spectroscopic target sample. There are clear differences: Quasars 
are in the most overdense environments at all scales, and Type II objects are in more over- 
dense environments than Type I objects for both higher luminosity AGN (i.e. quasars) and 
lower luminosity AGN. Within a scale of ~ 150 h^Q kpc, Type II quasars have an environ- 
ment 1.4 times more overdense than that of Type I quasars, albeit with large errors, while 
the Type II AGN have an environment 1.3 times more overdense than Type I AGN. At 
the same scale, the Type I quasars have environments more overdense than Type I AGN 
by a factor of 1.8, and Type II quasars have environments 1.9 times more overdense than 
Type II AGN. Moving out to the scale of ~ 1 Ji^q Mpc, the Type II quasars again have an 
environment 1.4 more overdense than the environment of Type I quasars, and the Type I 
quasars are in environments 1.4 times more overdense than Type I AGN. 

The differences between the target samples' environment overdensities could be an effect 
of AGN type, however, the intertwined effects of AGN luminosity and redshift will certainly 
play into these differences. In the next section, we explore how AGN type, redshift and 
luminosity influence the measured differences in AGN environments. 



4. DISCUSSION 

AGN luminosity and type should be key factors dictating variations in environment 
overdensity according to the hypothesis that information about the fueling mechanis ms of 



different AGN will be evidenced in the AGN environments (e.g.. iHopkins et al.ll2008l ). Ad- 
ditionally, the differences in overdensity between the target types seen in Figure [2] could be 
due to redshift effects, since it is obvious from Figure [1] that the redshift distributions of the 
samples are not the same. Therefore, in this section we investigate how the measurement of 
local environment varies with redshift, luminosity, and AGN type. We first isolate each of 
these variables to study its impact on measured overdensities while the other two variables 
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are held constant. Because the variables are not independent, we consider them together in 
Section H3J 

4.1. Redshift 

We first investigate only the effects of redshift. Figure [3] shows the mean cumulative 
overdensity as a function of redshift for the spectroscopic target samples, which provides 
marginal evidence for redshift evolution in the environment overdensity of Type I quasars. 
The magenta dashed line shows the linear weighted least-squares fit to the Type I quasar en- 
vironment overdensity data with redshift at different maximum radii; the fitting parameters 
for these lines is given in Table HJ The rightmost column of the table gives the % 2 proba- 
bility for each fit using the relevant degrees of freedom. The best linear fit for the smallest 
scale of R ^ 0.1 hjQ Mpc has a slope of —1.171, which suggests decreasing overdensity with 
increasing redshift, and for R ^ 0.25 Mpc, the overdensity increases with redshift: the 
best linear fit has a slope of 1.276. While these fits indicate a slight redshift dependence, 
we also try a zero-slope linear fit (which necessarily has one less degree of freedom than 
the best-fit line) and find that the zero-slope fit, i.e., no redshift dependence, is also a good 
fit to the data and in some cases is slightly more likely. Higher precision and higher red- 
shift measurements will be necessary to place strong constraints on the functional form of 
environment overdensity evolution with redshift. If the environment overdensity is indeed 
independent of redshift, this implies that the significant differences in environment seen in 
Figure [2] are caused primarily by luminosity and type effects, rather than the influence of 
redshift evolution. 

In the top panel of Figure HJ we show the evolution of the mean cumulative overdensity 
of photometric galaxies in the environments of Type I quasar, Type I AGN and Type II 
AGN samples. It is important to recall that we have placed the random points at the same 
redshift as the spectroscopic targets, and that we have imposed 5z cuts on the photometric 
galaxies (as described in Section [3D in order to minimize the effect of redshift evolution in 
the photometric galaxy sample. Therefore we can compare objects in different redshift bins. 
Figure H] demonstrates that higher redshift Type I quasars are in environments 1.24 times 
more overdense than the lower redshift quasars on scales < 500 kpc, while at larger scales, 
there appears to be little-to-no redshift evolution. However, there is scale-dependent redshift 
evolution evident on scales < 1.0 Mpc for the Type II AGN, shown in the lowest panel. 
The Type I AGN begin to exhibit noticeable redshift evolution at scales < 300 /if kpc, 
where the environments of lower redshift Type I AGN are 1.26 times less dense than those 
of the higher redshift Type I AGN. 
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We see therefore that there is some evidence for a change in local environment as a func- 
tion of redshift, all else being held constant. However, we have not yet taken AGN luminosity 
into account. Even in the same redshift range, selection effects due to the magnitude-limited 
samples may come into play, which we investigate in Sections 14.31 and 14.41 



4.2. Type 

In Figure [5l we identify three redshift ranges where there is overlap between our AGN 
samples and explore whether differences in type are reflected in the relative overdensity. The 
top panel shows the overdensity as a function of scale for both types of higher-luminosity 
AGN (i.e. quasars) in the range 0.3 ^ z ^ 0.6, and for both types of lower-luminosity 
AGN in two redshift ranges, 0.11 ^ z ^ 0.15 and 0.15 < z ^ 0.33. The dividing redshift 
value of z = 0.15 is chosen to roughly equalize the number of lower- luminosity AGN in each 
redshift range. The lower three panels show the ratio of Type II environment overdensity to 
Type I environment overdensity in the three redshift ranges. Again, we are able to compare 
objects in different redshift ranges because we have imposed 5z cuts on the photometric 
galaxies around both the spectroscopic targets and the random positions to which they are 
compared (see Section [3]) in order to account for any redshift evolution in the photometric 
galaxy sample and to minimize projection effects. 

Type II quasars have higher overdensity environments than Type I quasars with little 
scale dependence: at R « 1.0 h? Mpc, the overdensity of Type II quasar environments is a 
factor of 1.4 greater than the overdensity of Type I quasar environments, and at the smaller 
scale of R ~ 250 h^Q kpc, the Type II quasar environments have 1.2 times the overdensity 
of Type I quasar environments. However, the large errors prevent us from drawing strong 
conclusions. 

Type II environment overdensity is again consistently about a factor of 1.3 higher on 
all scales than the Type I environment overdensity for lower luminosity AGN in the redshift 
range 0.15 < z ^ 0.33. In the lower redshift range of 0.11 ^ z ^ 0.15, however, the Type 
II and Type I environment overdensities have a ratio consistent with unit y until scales R < 



200 hj n kpc, where the ratio increases to 1.3. This agrees with the work of iKoulouridis et al. 



( 120061 ). who found that Type II Seyfert galaxies were more likely to have a close neighbor 



than Type I Seyferts, using small samples of Seyfert 1 and 2 galaxies at very low redshifts. 

Because we see increased overdensity for Type II AGN compared to Type I AGN on 
small scales in overlapping redshift ranges, we can conclude that in Figure [2], the differences 
seen in environment overdensity between the AGN types are not primarily due to redshift 



- 11 - 



evolution. However, we have not ruled out the effects of AGN luminosity. The Type I and 
Type II AGN samples are selected from magnitude-limited spectroscopic galaxy samples, 
which will be dominated by intrinsically more luminous sources at higher redshift. The 
Type I AGN sample could be more affected by this magnitude limit, as the broad emission 
lines contribute more significantly to the overall flux in a given band and therefore the two 
AGN populations could have different average intrinsic luminosities. 



4.3. Luminosity 

U nlike Type I quasars, which are targeted largely based on their strong nuclear lumi- 



nosity (ISchneider et al.l 120071 ). the lower luminosity AGN we use were sel ected from objects 



classified as galaxies by the SDSS selection algorithms (lHao et al.l l2005al ) . The broad-band 
flux of these sources will be dominated by host galaxy starlight and/or flux from star for- 
mation, etc., which has little or no association with the nuclear luminosity. Therefore, we 
focus on the Type I quasar sample for our luminosity analysis, as it spans the entire redshift 
range we study, and with this long redshift baseline we are best able to disentangle redshift 
and luminosity effects on environment overdensity. 

In order to verify that the observed evidence for evolution of environment overdensity 
i s not due to the i ^ 19.1 (z < 3.0) limit imposed on Type I quasar selection in the SDSS 



(ISchneider et al.l 120071 ). we perform several tests in which we vary the apparent magnitude 
limit of the data. We consider two quasar samples limited to % ^ 18.9 and to the i ^ 19.1 
SDSS limit (see inset of Figure [1]). The two magnitude-limited samples were each then 
separated into two luminosity bins. We compared environment overdensity measurements 
of bright or dim quasars in each of the magnitude-limited samples and found no appreciable 
difference. Additionally, no difference was observed when different absolute magnitude values 
were used to define the bright and dim samples. In order to ensure that there is no difference 
between environments of quasars with i > 19.1, which were selected by the high-redshift 
targeting algorithm, and the rest of the apparent magnitude-selected sample, we performed 
similar tests comparing the environment overdensity of the entire quasar sample to that of 
the subset of quasars with i ^ 18.9 or i > 19.1. In all cases, there was no appreciable change 
in the observed overdensity. 

We compare the environment overdensities of Type I quasars in two luminosity bins to 
the other target samples without redshift cuts in Figure [6) The threshold value Mi = —23.2 
is chosen to give roughly equal numbers of Type I quasars in each luminosity bin: there 
are 2,190 (2,044) quasars with -27.5 ^ M { ^ -23.2 (-23.2 < M { ^ -22.0). The average 
magnitude of the brighter (fainter) bin is Mj = —23.83 (Mj = —22.70). 
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Type II quasars and the brighter Type I quasars are located in similarly overdense 
environments consistently at all scales, while the dimmer Type I quasars are located in 
environments slightly less overdense than the Type II quasars. At a scale R ~ 500 Hjq kpc, 
the cumulative overdensity of Type II quasar environment is 1.06 times that of the brighter 
Type I quasars, but 1.3 times as the dimmer Type I quasars. At the scale of R ~ 1.0 Mpc, 
Type II quasars have environment overdensities 1.2 times the environment overdensity of 
brighter Type I quasars but 1.5 times that of dimmer Type I quasars. Again we note that 
the large error bars nearly overlap with unity and prevent strong conclusions. 

The more luminous Type I quasars are located in environments more overdense than 
Type I AGN, while there is less difference in the overdensities of dimmer Type I quasars 
and Type I AGN. The environment overdensity ratio increases with decreasing scale for both 
brighter and dimmer Type I quasars. At a scale R ~ 500 h? kpc, brighter Type I quasar en- 
vironments have an overdensity 1.6 times the overdensity of Type I AGN environments with 
significance ~ 3a, and dimmer Type I quasar environments have an overdensity 1.3 times 
the overdensity of Type I AGN environments with significance m 2a. At R fa 150 h^Q kpc, 
the environments of brighter Type I quasars are 2.1 times as overdense (2.4a), and the envi- 
ronments of dimmer Type I quasars are 1.6 times as overdense as the environments of Type 
I AGN (2a). 

The ratio of Type I quasars to Type II AGN increases for both bright and dim quasars 
with decreasing scale, but less dramatically as the ratio to Type I AGN. The ratio between 
dimmer Type I quasars and Type II AGN is approximately consistent with unity for scales 
150 Jiyq kpc < R ^ 2.0 Mpc; the ratio between brighter Type I quasars and Type II 
AGN is 1.3 (> 2a) for scales R > 500 h^Q kpc. On smaller scales, both ratios increase. 
At scales R ~ 150 kpc, the ratio of brighter Type I quasars to Type II AGN is 1.6 
(^ 2a), and the ratio of dimmer Type I quasars to Type II AGN is 1.3 (> la). This scale 
dependency could be evidence for the merger origin of quasars, since one would expect to 
see a higher d ensity of environmen t galaxies at small scales where merger events are likely 



to take place (IHopkins et al.ll2008l ). 



Evidence that dimmer quasars and lower-luminosity AGN are located in environments 
with similar overdensity might suggest that dimmer quasars could be a transition popula- 
tion between low-luminosity AGN (likely fueled in dry mergers, close encounters, or secular 
processes) and high- luminosity AGN (likely fueled in major mergers). Rather than dis- 
parate populations of merger-fueled and secularly fueled AGN, there may be a continuum of 
galaxy interactions from major mergers to close encounters or harassment that cause AGN 
luminosity differences. Alternatively, a mix of mergers and secula r processes could drive 



the AGN population near the quasar-Seyfert divide (Mj « —22.5; iHao et al.ll2005bl ). We 
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have compared the AGN samples without redshift cuts, but we note that in Section 14.11 we 
demonstrated that evolution of quasar environments with redshift is negligible. 

The significant difference in the environments of bright Type I quasars and the environ- 
ments of both Type I and Type II AGN could imply that these pop ulations have d i fferen t 
fueling mechanisms. This result is consistent with results presented by iLi et al.l (120061 . 120081 ) , 
who find that there is only a weak link between nearby neighbors of narrow-line AGN and 
their nuclear activity. 



4.4. Luminosity, Redshift, and Type 



Finally, we combine our analysis of type, redshift and luminosity effects on environment 
overdensity in Figures [7J [HJ and [9j Our 8z cuts on the photometric galaxies around the 
spectroscopic targets as well as around the random positions to which they are compared 
(as described in Section l3TTj) allow us to make meaningful comparisons of objects in different 
redshift ranges. In Figure 0, Type II quasars are compared to Type I quasars in the redshift 
range 0.3 ^ z ^ 0.6. We divide the Type I quasars into bright (2,001; Mj = —23.87) and 
dim (1,915; M, = —22.75, about 2.8 times fainter) samples of roughly equal numbers at 
Mi = -23.25. 

Comparing the lower panel of Figure [7] to the top ratio panel of Figure [5] shows the 
dramatic part luminosity plays compared to evolution alone. The environment of Type II 
quasars is similar to the signature of brighter Type I quasars for all scales we measure. The 
similarity of environments at small scales suggests that the differences observed between 
brighter Type I quasars and Type II quasars are due to a non-environmentally driven mech- 
anism such as orientation or internal structure effects. This in turn implies that Type II 
quasars are not a different cosmological population from these brighter Type I quasars. 

However, the Type II quasar environments are consistently more overdense, albeit with 
large errors, than those of dimmer Type I quasars on scales R ^ 2.0 /ifg 1 Mpc that we 
measure. The different characteristics of the environments of the dimmer Type I quasar 
population from the Type II quasar population are most likely due to intrinsic luminosity 
differences rather than redshift differences. We see consistent overdensity ratios on all scales 
and do not see small scale effects, therefore we conclude that the difference in environment 
overdensity between the brighter and dimmer quasars is primarily du e to mass effects. More 
luminous AGN are ex pected to have higher mass black holes (e.g., iMagorrian et al.l 11998 



le.g. 



Marconi fc Hunt 2003). which are in turn correlated with more massive dark matter halos 



Ferrarese fc Merrittll2000l ; iGebhardt et all2000l ; lTremaine et al.ll2002l ). Selection effects 
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in the magnitude-limited photometric galaxy sample could also play into the difference in 
overdensity between the brighter and dimmer Type I quasars. The redshift distribution 
of Type I quasars in the brighter (Mj ^ —23.25) bin is slightly different from that of the 
dimmer (M, > —23.25) bin even over the redshift range of 0.3 ^ z ^ 0.6 (the brighter quasars 
have a mean redshift of 0.51, and the dimmer quasars have a mean redshift of 0.45). The 
galaxies seen in the environments of brighter (higher redshift) quasars will tend themselves 
to be brighter, and consequently more massive, and therefore cluster mo re strongly than 



the dimmer environment galaxies (IMaddox et al.lll990l ; IZehavi et al.ll2002l ). This, however, 
should not be a major effect. 

Figure [8] compares the environments of Type I quasars in two luminosity bins to the 
environments of Type I and II AGN in the redshift range 0.15 < z ^ 0.33. We use Mj = 
—22.65 as the threshold value for brighter and dimmer quasars in this lower redshift range to 
equalize the number in each luminosity bin. The 244 brighter (239 dimmer) Type I quasars 
have a mean magnitude of Mj = —22.33 (Mj = —23.30, about 2.4 times fainter than the 
brighter sample). We note that these Type I quasar samples are about a factor of seven times 
smaller than the Type I quasar samples in the higher redshift range, thus the measurements 
(and resulting interpretation) will be less precise. 

The top panel of Figure [8] shows that at for scales R > 300 h^Q kpc, the environments 
of dimmer Type I quasars are more overdense than those of brighter Type I quasars. It 
appears that the situation has been reversed from Figure where the environments of 
dimmer quasars were less overdense than the environments of brighter quasars. However, 
the range of luminosity at this lower redshift range of 0.15 < z ^ 0.33 is much smaller than 
for the higher range of 0.3 ^ z ^ 0.6. The overall absolute magnitude distribution of these 
Type I quasars is skewed toward the faint end, thus the divi ding value Mj = — 22.65 is very 



close to the quasar-Seyfert divide of Mj ps —22.5 as defined by lHao et al.l (l2005bl ). Significant 
variation in overdensities was seen when different magnitude cuts were imposed, with the 
dimmer quasars consistently having higher overdensities by varying margins. The dramatic 
sensitivity of results on the bright/dim dividing value emphasize that for low luminosites 
and redshifts, broad-band absolute magnitudes are a poor proxy for AGN luminosity. The 
measured flux is more likely to be affected by galaxy starlight, star formation, etc. at 
this faint end. Therefore, any attempt to use broadband magnitudes to correlate nuclear 
luminosity with environment will be skewed. 

With these caveats in mind, we compare the Type I quasars to lower-luminosity Type I 
and Type II AGN in the lower two panels of the figure. Dimmer Type I quasar environments 
have overdensities greater than the Type I AGN, but the environments of brighter Type I 
quasars and the Type I AGN have about the same amplitude on all scales. The lower ratio 
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panel shows the ratio of bright and dim Type I quasars to Type II AGN. The brigher quasars 
have environments with slightly lower overdensity than Type II AGN; the environments of 
dimmer Type I quasars are only slightly more overdense than the environments of Type II 
AGN, but are consistent within the error bars. 

In Figure [9] we focus on Type I quasars alone to investigate the evolution of the environ- 
ment overdensity of brighter and dimmer objects. We have chosen two luminosity intervals 
of one magnitude in width and compare the environment overdensities of brighter to dimmer 
objects in three redshift intervals. AMi corresponds to the dimmer luminosity interval of 
—23.0 < Mi ^ —22.0 and contains 1,556 Type I quasars, and the brighter luminosity interval 
AM2 is —24.0 < Mi ^ —23.0, containing 2,036 Type I quasars. Table [2] gives the number of 
quasars as well as the mean magnitude in each redshift and magnitude bin. 

In the two lower redshift bins 0.15 < z ^ 0.3 and 0.3 < z ^ 0.45, there is little 
difference in the environment overdensity of brighter and dimmer quasars with little-to-no 
scale dependence. However, in the highest redshift interval of 0.45 < z ^ 0.6, brighter 
quasars are shown to be located in slightly more overdense environments than the dimmer 
quasars. At scale of R ~ 1.0 h^Q Mpc, the brighter quasars are located in environments with 
overdensity 1.5 times that of the dimmer quasars. The brighter quasars have environments 
with overdensity 1.4 times the overdensity of dimmer quasar environments at a scale of 
R w 250 hjQ kpc, and then the ratio begins to drop toward unity at the innermost scales. 
However, the large errors are nearly consistent with unity on all scales we measure. 

It appears, therefore, that there is again slight evidence for some redshift evolution of 
Type I quasar environments, but it is mainly manifested at the highest redshift range. This 
emphasizes the need for additional studies of the environments around higher redshift Type 
I quasars. We caution that the increased overdensity at higher redshift may be influenced by 
the fact that there are nearly three times the number of brighter quasars as dimmer quasars 
in the interval 0.45 < z ^ 0.6 (see Table [2]). In the range 0.3 < z ^ 0.45, the number of 
bright quasars is closer to the number of dim quasars, while in the range 0.15 < z ^ 0.3, 
the dim quasars outnumber the bright quasars by more than a factor of two. The change 
in overdensity ratio will also be affected by the change in mean luminosity of the dimmer 
quasar sample with increasing redshift. While the bright quasar luminosity changes only by 
0.07 magnitudes, the mean dim quasar luminosity changes by 0.17 magnitudes. Therefore 
we cannot draw strong conclusions, but reiterate the need for higher precision and higher 
redshift measurements of quasar environments. 
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5. CONCLUSIONS 



Our work sheds new light on the nature of AGN environments and their relationship to 
the type, luminosity, and redshift of the AGN itself. We have used larger samples of AGN 
targets and imposed a photometric redshift cut on the nearby photometric galaxies in order 
to minimize projection effects and to account for redshift evolution of the photometric galaxy 
sample. By using photometric redshift cuts, we obtain more realistic overdensity estimates 
and errors for the local environments of AGN of various luminosities and types. 

There are two main pictures through w hich observed differen ces in AGN can be inter- 
preted. In the merger models presented by iHopkins et al.l ( 120061 ) . differing fueling mecha- 
nis ms trigger AGN a cti vity at differe nt luminosities. Unified models such as those presented 
by lAntonuccil (119931 ) or lElvid (120001 ) ascribe observed differences in AGN to viewing angle 
or structure. These pictures are not mutually exclusive; in fact, we show that the interde- 
pendency between the variables of type, luminosity and redshift play into the subtleties of 
distinguishing the pictures based on environment overdensity measurements. 

From our analysis, we can draw the following conclusions: 



• Type II quasars are shown to have similar environments as brighter Type I quasars in 
the same redshift range on all scales that we study, which suggests the observational 
differences in Type I and Type II quasars are driven by orientation and/or structure 
and not by cosmological evolution. 

• Lower-luminosity Type II AGN have environments consistently more overdense than 
lower-luminosity Type I AGN by a factor of 1.3 at all scales in the redshift range 
0.15 < z ^ 0.33, but there is scale dependence in the environment overdensity ratio of 
Type II to Type I AGN for redshifts 0.11<2<0.15. 

• There is marginal evidence for redshift evolution of Type I quasar environments on 
all scales, especially for 0.45 ^ z ^ 0.6, not noted in previous studies. However, this 
evolution is not the primary explanation for the environment overdensity differences 
seen between Type I quasars and Type II quasars, and between Type I AGN and Type 
II AGN. In order to place strong constraints on the functional form of this redshift 
evolution, it is necessary to acquire higher precision measurements and higher redshift 
measurements. 



In our current work, we have explored the effects of redshift, type and luminosity on the 
measurement of AGN environments in order to better understand the evolutionary nature 
of AGN and their general relationship with galaxies. There are several ways in which our 
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existing work can be improved. First, when we measured the redshift effects of the AGN 
environments, we were confronted by the selection effects that result from our use of a 
magnitude-limited photometric galaxy sample. The standard way to overcome this effect 
is to employ volume-limited samples of galaxi es, which can be done by using photometric 
redshifts to derive absolute magnitudes (e.g., iRoss et all 120071 ). This approach provides a 
more accurate technique to compare AGN environments across different redshift ranges. 
Second, as we have discussed in Section 14.41 broad-band absolute magnitudes are not the 
best way to measure the nuclear luminosity of AGN. Accretion luminosity, which can be more 
easily compared to the predictions of quasar formation models, can be better quantified by 
using emission line luminosities, such as 0[III]. Third, our current analysis simply studies 
the local environment of AGN by measuring the overdensity of all galaxies near the AGN. 
This approach, while useful, ignores a great deal of information. The physical characteristics 
of the galaxies in the environment pr ovide additional yaluable insight into the re lationship 
between AGN and galaxies (see, e.g. iBarr et al.ll2003l ; IColdwell et al.ll2003l 120061 ) . Finally, 
by exploring the relationship between AGN environments and intrinsic physical properties 
of the AGN, such as black hole mass and multiwavelength emission, we can further refine 
our understanding of AGN fueling mechanisms. 
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Table 1. Linear least-squared fit parameters for quasar environment data in Figure 



Rmax slope intercept \ 2 P(xV) 



2.0 hn Mpc 0.119 -0.880 6.84 0.3363 
0.0 -0.833 8.10 0.3241 
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Table 2. Details for data used in Figure [9j 



Redshift Range 


-24.0 < Mi < -23.0 


-23.0 < Mi < -22.0 


0.15 < z ^ 0.3 


85, -23.39 


211, -22.46 


0.3 < z ^ 0.45 


502, -23.39 


788, -22.57 


0.45 < z ^ 0.6 


1445, -23.46 


552, -22.63 
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Fig. 1. — Redshift distribution of spectroscopic targets and photometric galaxies. The 
photometric galaxy redshift distribution is shown by the solid black curve. The Type I 
(—27.5 ^ Mi ^ —22.0) point-source quasar distribution (4,234 objects) is shown with the 
black long-dashed histogram, and the Type II quasar distribution (160 objects) with the red 
dotted histogram. The distribution of lower-luminosity Type I AGN (1,745) is shown with 
the green short dashed histogram, and the Type II AGN distribution (classified by Kewley's 
criteria; 3,408) is shown with the blue dot-dashed histogram. The inset shows absolute i- 
magnitude of Type I quasars v s. redshift. Evidence f or the i ^ 19.1 (z < 3.0) limit used in 
the quasar selection algorithm (ISchneider et al.l 120071 ) is seen as a thicker diagonal band. 
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Fig. 2. — Mean cumulative overdensity of photometric galaxies as a function of comoving 
scale around spectroscopic targets. Solid black squares represent Type I quasars, open red 
starred points represent Type II quasars, solid blue triangles represent Type II AGN, and 
solid green hexagons represent Type I AGN. Points have been offset slightly for clarity. 
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Fig. 3. — Mean cumulative overdensity vs. redshift for spectroscopic targets. Symbols 
correspond to those used in Figure [2j The magenta dashed lines are linear weighted least- 
squares fits for the Type I quasar sample; the parameters for these lines are given in Table [TJ 
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Fig. 4. — Upper panel: Mean cumulative overdensity of photometric galaxies around quasars 
and lower-luminosity AGN as a function of scale and redshift. Points have been offset slightly 
for clarity. Top lower panel: Ratio of environment overdensity of higher-redshift Type I 
quasars to that of lower-redshift Type I quasars. Middle lower panel: Ratio of environment 
overdensity of higher-redshift Type I AGN to that of lower-redshift Type I AGN. Bottom 
lower panel: Ratio of environment overdensity of higher-redshift Type II AGN to that of 
lower-redshift Type II AGN. 
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Fig. 5. — Upper panel: Mean cumulative overdensity of photometric galaxies around quasars 
and lower-luminosity AGN as a function of scale and redshift. Points have been offset slightly 
for clarity. Top lower panel: Ratio of environment overdensity of Type II quasars to that of 
Type I quasars in the redshift range 0.3 ^ z ^ 0.6. Middle lower panel: Ratio of environment 
overdensity of Type II AGN to that of Type I AGN in the redshift range 0.15 < z ^ 0.33. 
Bottom lower panel: Ratio of environment overdensity of Type II AGN to that of Type I 
AGN in the redshift range 0.11 < z ^ 0.15. 
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Fig. 6. — Top panel: Mean cumulative overdensity of photometric galaxies around Type I 
point-source quasars split by luminosity; low luminosity AGN; and Type II quasars. Bright 
quasars have absolute magnitude —27.5 ^ M» ^ —23.2 and dim quasars have absolute 
magnitude —23.2 < Mi ^ —22.0. Top lower panel: ratio of Type II quasar environment 
overdensity to bright (solid points) and dim (open points) Type I quasar environment over- 
densities. Middle lower panel: ratio of bright (solid points) and dim (open points) Type I 
quasar environment overdensities to Type I AGN overdensity. Bottom lower panel: ratio 
of bright (solid points) and dim (open points) Type I quasar environment overdensities to 
Type II AGN. No redshift limits have been imposed. 
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Fig. 7. — Upper panel: Mean cumulative overdensity of photometric galaxies around Type 
II quasars Type I quasars in the redshift range 0.3 ^ z ^ 0.6. The Type I quasars in this 
redshift range have been divided at Mj = —23.25 so that the luminosity bins contain approx- 
imately equal numbers of Type I quasars. Lower panel: Ratio of environment overdensities 
of Type II quasars to brighter and dimmer Type I quasars in this redshift range. Points in 
both panels have been offset slightly for clarity. 
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Fig. 8. — Upper panel: Mean cumulative overdensity of photometric galaxies around Type 
I quasars and Type I and II AGN in the redshift range 0.15 < z ^ 0.33. The Type I 
quasars in this redshift range have been divided at Mj = —22.65 so that the luminosity bins 
contain approximately equal numbers of Type I quasars. Lower panels: Ratio of environment 
overdensities of brighter and dimmer Type I quasars to Type I and Type II AGN in this 
redshift range. Points in both panels have been offset slightly for clarity. 
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Fig. 9. — Mean cumulative overdensity of photometric galaxies around Type I quasars with 
redshift for two luminosity bins, where AM 1 : -23.0 < Af< < -22.0 and AM 2 : -24.0 < 
Mi ^ —23.0. Lower panels: Ratio of AM 2 to AMi quasar environment overdensities in the 
three redshift ranges. 



